Enhanced Josephson tunneling between high temperature superconductors through a normal 
pseudogap underdoped cuprate with a finite energy cooperon 
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The Josephson coupling between optimally cuprate superconductors separated by a spacer with a finite en- 
ergy cooperon excitation which contributes to the Josephson coupling strength, is examined. For an underdoped 
cuprate barrier in its normal state, the YRZ model gives a good description of the temperature dependent en- 
hanced Josephson coupling. A detailed examination of origin of the enhancement shows a significant contribu- 
tion from the cooperon excitation which is comparable to that from nodal quasiparticles. 
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There have been a number of reports of enhanced, even 
giant, penetration lengths for d-wave superconductivity into 
nonsuperconducting phases of underdoped cuprates 1-4 . Bo- 
zo vie et al. reported giant proximity effect in uniform trilayer 
junctions of high temperature superconductors at temperature 
T above the superconducting transition temperature T c of the 
spacer underdoped cuprate layer in the pseudogap phase 2 . The 
enhanced penetration length observed in this case has been 
explained by a model with strong phase fluctuations due to a 
reduced phase stiffness at underdoping 5 . 

Recently it has been proposed that the opening of the 
pseudogap in the antinodal regions of &-space in underdoped 
cuprates, is accompanied by the appearance of a finite energy 
bound hole pair excitation, also known as a cooperon 6 ^. In 
this letter we examine the influence of such cooperon excita- 
tions on the penetration length of superconducting order into 
nonsuperconducting cuprate materials. A cooperon creates a 
low energy pole in the Cooper pair correlation function lead- 
ing to enhanced penetration of superconducting order. We 
show that the cooperon contribution to the enhancement of 
the Josephson tunneling in thick layer is significant. 

We are motivated by the recent ARPES 8-11 (angle resolved 
photoemission spectroscopy) and other experiments 12 ~i2i22, 
which support rather a 2-gap scenario with an insulating pseu- 
dogap for underdoped cuprates. This leads us to examine 
the superconducting penetration at T > T c , using the phe- 
nomenological 2-gap YRZ theory put forward by Yang et al. 
for the pseudogap phase. The YRZ model has had consid- 
erable success in describing many anomalous properties of 
the pseudogap phase. For a recent review see Rice,Yang and 
Zhang 21 . As discussed below, this model has a low energy 
cooperon pole at temperatures T > T c which enhances the 
penetration. 

An analysis of the single particle propagator in a 2- 
dimensional array of 2-leg Hubbard ladders was an important 
input in the formulation of the YRZ ansatz. A recent exten- 
sion by Konik, Rice and Tsvelik to an array of 4-leg Hubbard 
ladders with 4 bands crossing the Fermi surface. They showed 
that at low doping there is an insulating gap with a cooperon 
resonance in the outer band pair which coexists with metal- 



lic behavior in the inner band pair. The virtual exchange of 
the cooperon resonance with support on the outer bands intro- 
duces a pairing mechanism in the inner band Fermi surface 6 . 

A decade ago Nozieres and Pistolesi 22 studied the transition 
for fermions with an attractive interaction between a super- 
conducting and a semiconducting state within BCS mean field 
theory as the one particle band gap is increased. They found 
that at a critical value of the band gap, which is simply related 
to the superconducting gap of the starting metallic state, the 
pairing amplitude and phase stiffness dropped continuously to 
zero while the single particle gap remained finite. Recently 
Rice et al— showed that the approach to this quantum critical 
point(QCP) from the semiconducting state was characterized 
by a softening of the cooperon mode, dropping to zero en- 
ergy at the QCP 21 . Support for the existence of cooperons in 
strongly underdoped cuprates also comes from exact diago- 
nalization results of the strong coupling t-J models on small 
clusters. When extrapolated to an infinite lattice, these show a 
low lying cooperon resonance in the antinodal regions of the 
Brillouin zone 6 . 

We use the phenomenological YRZ theory to describe the 
pseudogap phase with input parameters chosen to have the 
same values as in the original work 8 to model the experi- 
ments by Bozo vie et. al. on the enhanced superconducting 
penetration 2 into underdoped LSCO cuprates at temperatures 
T > T c . In the YRZ model, the coherent part of the Green's 
function is given by, 
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The self-energy is given by £#(k, co) = \A R (k)\ 2 /[co + e(k)] 
with 6k = -2t(cosk x + cos k y ), £k = £k ~ 4?' cos k x cos k y - 
2t"(cos2k x + cos 2k y ) - /u p , and the RVB gap A R (k) = 
A (x)(cos k x - cos £3;). g t = 2x/(l + x) is the doping-dependent 
renormalization factor 23 , t, t\ t" are the renormalized hop- 
ping integrals, and fi p is an effective chemical potential deter- 
mined by Luttinger sum rule that the area enclosed by the four 
Fermi pockets equals to the doped hole density 8 ' 24 . The modi- 
fied pairing form for the self-energy leads to two quasiparticle 
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bands with energy dispersion and the corresponding spectral 
weight, 



£k = (fik - ek)/2 ± ^(fik + ^k) 2 /4 + A 2 (k) 
K = [l+(fik + ^)/(2Ej-fik + ^)]/2 



(2) 
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Choosing a d- wave attractive interaction Vi(k||,q||) = 
ViJ(k||)d(qu) with V t < 0, then the self-consistent 
superconducting-gap equation is given by, 



Ai = |Vi|Im2 J ^J(k|,)/^)[F^k||,^)-F?;-(k||,^)] (6) 



with /? = l/(k B T) the reciprocal temperature, /(jc) = l/(e x + 
1) the Fermi function, F r/fit the retarded/advanced anomalous 
Green's function. According to Eq.©, gap equation Eq.© 
can be simplified as, 



Ai = T Z Z ^(kn)/0S^ k|| )< k|| (/)[v- k|| (/)]* (7) 



FIG. 1 : Schematic illustration of the setup. The leads are two optimal 
doped LSCO material, while the center is an underdoped LSCO with 
doping x = 0.1. The interface is parallel with the CuO^ plane. 

We introduce superconductivity by adding a pairing poten- 
tial with d-wave symmetry. The magnitude of the pairing po- 
tential is chosen to reproduce mean-field values for T c in the 
Bozovic experiments, i.e. T c = 25K (doping x = 0.1) in the 
underdoped junction and T c = 45 K (doping x f = 0.15) in the 
optimum-doped leads. The tunneling setup is represented by 
a set of layers, typically 10 for each of the leads and 10 for the 
junction, as shown in fig. \T\ Next we obtain a set of coupled 
Bogoliubov-de Gennes (BdG) equations for the layers, 
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where the BdG Hamiltonian Hf. is given by, 



£?(k|l)tfy - t±M , Af(k||)5y 

A*(k„)^, -£f(k||>5y + f^(k||) 



(4) 



with i,j the layer index, ky the in-plane momentum, Ef(k\\) 
the YRZ quasiparticle dispersion given by Eq.© which de- 
pends on the doping value of of layer-/. £ ± ^-(ky) = t c (cosk x - 
cos k y ) 2 / '4 is the standard hopping between nearest-neighbor 
layers 25 (//), and A/(ky) = A;d(ky) is the superconducting gap 
of layer-/ with the d- wave factor d(ky) = cosk x - cosk y . 8® k 

and [w^ k|| , v^ k| ] r are the eigenvalues and the corresponding 
orthonormal eigenvectors of BdG equation. 

The Nambu Green's function for this multi-layer system 
can be constructed by the wave-functions as 2 ^, 
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The Josephson current between nearest-neighbor layers is 
jiven by 



hi = (2e/») J] KiM^i^hM + c lu c hM> (8) 



This can be evaluated as 27 , 



(9) 

Iij satisfies the continuity equation 28,29 between the leads, i.e., 
independent of layer index. Fig.© gives the result of the crit- 
ical Josephson current as a function of temperature obtained 
from YRZ model. We choose the thickness of underdoped 
LCO d ~ 10 layers, the perpendicular hopping t c = 0.05^o 
with to ~ 300meV&22, the doping x = 0.1, the magnitude 
of d-wave attractive interaction V = -0.54^, which gives 
a mean-field critical temperature around 25^ as in the Bo- 
zovic experiments 2 . We find a decay length of the Josephson 
current around A ~ 10(unit cell) at a temperature T = 30 K 
which is higher than the critical temperature of the under- 
doped spacer T c « 25 K. This longer range penetration at 
temperatures above the superconducting transition of the un- 
derdoped spacer is due to an enhanced proximity effect. In 
contrast, for an insulating spacer without cooperon states (in- 
duced by attractive interaction), a single layer is enough to kill 
the Josephson coupling in our formulism. These results agree 
with the experiments 2,31 on the stoichiometric LCO quite well. 

The cooperon excitation in underdoped LSCO at tempera- 
tures T > T c with total momentum q z , appears as the pole at 
an energy in the pair propagator which is defined in real 
space as, 



3 




FIG. 2: The critical Josephson current as a function of temperature. 
The parameters for the underdoped spacer are doping x = 0.1, the 
magnitude of perpendicular hopping t c = 0.05 to, the magnitude of 
d- wave interaction V = -0.54^. The black curve is for the exper- 
imental data. The inset gives the cooperon energy as a function of 
temperature obtained from Eq. (ITTi with the same parameters. 
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FIG. 3: The ladder-diagram expansion of the pair propagator and 
its relation to Josephson tunneling. Upper figure: Feynman diagram 
of Josephson tunneling; Lower figure: ladder-diagram expansion of 
the pair propagator. The solid line denotes the one-particle Green's 
function, and the wavy line denotes the effective d-wave interaction, 
V' qk = V'd(k)d(q). Note that V'= V/Z, with V the bare d-wave in- 
teraction and Z the spectral weight. See also Schachinger and Car- 
botta^ 



Its Fourier transform in momentum space can be estimated 
by summing over the particle ladder diagrams 33,34 See Fig. 
©, which gives the result, 



G 2 (0, io) n ) 



(10) 



when the total momentum q z = and the perpendicular hop- 
ping t c is much smaller than the in-plane hopping to. Here 
the factor X° (ia) n ) = tanh(J3E%J2)/[2E^ + irj - ioj n ] with 
and denotes the dispersion and spectral weight of YRZ 
dispersion given in Eq.©. rj = 0.01 to is a small imaginary 
part, representing the decay of the cooperon state 6 . The zero- 
momentum cooperon energy Qo is determined by the real part 
of the pole of pair propagator given by Eq. (Hoi i.e., 
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The insert in Fig.© shows the decrease (in electron nota- 
tion^ Do from zero as T increase above T c , in line with the 
behavior seen in numerical simulations 7 of a related model. 

The contribution of the cooperon excitation can be deduced 
from the distribution of the Josephson current in momentum 
space of the underdoped sample (doping x « 0.1). In the 
antinodal region of the Brillouin zone where cooperon states 
play an important role 6 , the Josephson current is enhanced by 
cooperon states. On the other hand, Josephson current can 
also be contributed from the quasi-particle excitation around 
the Fermi pockets centered at nodal region of Brillouin zone 8 . 
The right panel of Fig.© gives the distribution of the Joseph- 
son current in Brillouin zone, which peaks at the "banana tips" 
of the Fermi pocket. Then we calculate the contributions from 
nodal region, the region between the two white lines in the 
right panel of Fig 01 and the antinodal region, the region out- 
side the white lines. Though the contribution at a given k in 
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FIG. 4: Left panel: The contribution to the Josephson current from 
nodal and antinodal region, respectively. Right panel: The distribu- 
tion of Josephson current in Brillouin zone. The choice of parame- 
ters are same as the ones in Fig|2] The region between the two white 
lines are nodal region, while the others are antinodal region. Tough 
the Josephson current is strongest at the"banana tips" of Fermi pock- 
ets, the total contribution from nodal region is comparable with the 
one from broadly spreading current around the antinodal region. 

antinodal region is smaller than the "banana tips", the rela- 
tively wide area in the antinodal region makes its total con- 
tribution to the tunneling comparable to the contribution from 
the nodal region including the "banana tips" as shown in the 
left panel of fig.H Note that in LSCO the neighboring Cu0 2 
planes along the c-axis are displaced by half a lattice param- 
eter which reduces the weighting of the antinodal regions 20 
in t ± jj(k\\) relative to the standard interlay er hopping matrix 
element and the contribution of the cooperon excitation. 

There are also a few experiments on 
YBCOI Y x Pri- x Ba 2 Cu 3 7 / YBCO superstructure, which 
also shows a very large proximity length at x < 0.55 1,4 . The 
spacer material YPBCO is semiconducting at x = and only 
becomes superconducting with increasing Y- substitution at 



x > x c ( = 0.5). Recently Wojek and collaborators 4 used 
low-energy muon-spin rotation(LE-yuSR) to measure the 
pentration of the superconducting order parameter into a 
stoichiometric insulating PBCO spacer. They reported a 
substantial penetration of the superconductivity with a large 
spacer width of 45nm. Both the contributions from nodal 
quasiparticle and cooperon are weak in this case, and the 
strong proximity effect needs a different explanation. 

The phenomenological YRZ theory for the pseudogap 
phase at underdoping interprets the antinodal energy gap as 
a precursor to the Mott insulator at zero doping. The analogy 
between the 2-dimensional and finite ladder systems has led 
to proposals that a cooperon resonance appears as the antin- 
odal energy gap opens up. The aim of this paper is to test this 
proposal by examining the Josephson coupling between opti- 
mally (or overdoped) cuprate superconductors separated by a 
barrier of an underdoped cuprate in its pseudogap state with 
a finite energy cooperon excitation which contributes to the 
Josephson coupling strength. Here the calculations based on 
the YRZ model give a good description of the temperature de- 
pendent enhanced Josephson tunneling. A detailed examina- 
tion of origin of the enhancement shows that the contribution 
from the cooperon excitation is significant and is comparable 
with that from the nodal quasi-particles. 
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